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At the Fcrmilab Tevatron, where pp collisions are pro- 
duced at y/s = 1.96 TeV, W + (W~) bosons are created 
primarily by the interaction of u (d) quarks from the pro- 
ton and d (u) quarks from the anti-proton. Since u quarks 
carry, on average, a higher fraction of the proton's mo- 
mentum than d quarks [HQ > ^ ne ^ + tends to be boosted 
along the proton beam direction and the W~ tends to be 
boosted along the anti-proton direction. The difference 
between the W + and W~ rapidity distributions results 
in a charge asymmetry 



A(y w ) 



da + 1 dyw — da /dyw 
da + 1 dyw + da^ /dyw 



(1) 



where yw is the W boson rapidity Q and da^ j dyw is the 
differential cross section for W + or W~ boson produc- 
tion. The parton distribution functions (PDFs) describ- 
ing the internal structure of the proton are constrained 
by measuring A{yw) 01- 

Previous measurements [D, S 0, S] of the W charge 
asymmetry at the Tevatron were made as a function of 
the pseudorapidity r\ 0| of the leptons from decays of 
W — ► lui (I = e, /it) since the W decay involves a neu- 
trino whose longitudinal momentum is not determined 
experimentally. However, the lepton charge asymmetry 
is a convolution of the W production charge asymme- 
try and the V — A asymmetry from W decays. These 
two asymmetries tend to cancel at large pseudorapiditics 
(M ^ 2.0), and the convolution weakens and complicates 
the constraint on the proton PDFs. 

In the measurement presented in this Letter, the com- 
plication is resolved by using additional information in 
the lepton transverse energy (E T ) and the missing trans- 
verse energy ($ T ) [3 on an event-by-event basis to mea- 
sure the asymmetry as a function of the \yw\ instead 
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of the lepton \r/\. This new analysis technique [9J] gives 
the first direct measurement of the W production charge 
asymmetry using W — ► ev decays. We use data from 
1 fb _1 of integrated luminosity collected by the CDF II 
detector. The region of acceptance is \yw\ < 3.0, giv- 
ing the new measurement an ability to improve proton 
PDFs determinations for 0.002 < x < 0.8, where x is the 
fraction of the proton momentum carried by u- or d-typc 
quarks. This analysis is described in detail in (loj . 

The CDF II detector is described in detail elsewhere 
What follows is a brief description of the detector 
components needed to identify W — > ev events, which 
are characterized by large missing transverse energy (j£ T ) 
and a track in the central drift chamber (COT) [12j 
or in the silicon tracking system (SVX) [l3|, [l4| that 
points to a cluster of energy in the electromagnetic (EM) 
calorimeters [ljl [lj|- The SVX provides precise track 
measurements from eight radial layers of microstrip sen- 
sors. The COT provides additional tracking information 
from 96 layers of wires. Tracks are measured inside a 1.4 
T solenoidal magnetic field that allows electron charge 
determination from the curvature of the track. The 
COT allows track reconstruction in the range \rj\ < 1.6, 
while the SVX extends the capability up to \r)\ ~ 2.8. 
Outside the tracking system, EM and hadronic (HAD) 
calorimeters measure the energies of showering particles. 
The calorimeters are divided into two types: a central 
calorimeter with a fiducial region covering \rj\ < 1.1, and 
a forward calorimeter covering 1.2 < \rj\ < 3.5. 

We use two types of W — > ev events, classified by the 
calorimeter section in which the electron is detected. The 
data are initially selected by an on-line event selection 
(trigger) system. The trigger for the central electrons 
requires an EM energy cluster with E T > 18 GeV and a 
matching track with p T > 9 GeV. The forward trigger, 
designed specifically for W candidates, requires an EM 
energy cluster with E T > 20 GeV and £ T > 15 GeV. 

For central electrons, we require off-line event selection 
including an isolated energy cluster in the region |ry| < 1.1 
with E T > 25 GeV and Iso(0.4) < 4.0 GeV. The isolation 
Iso(0.4) is defined as the calorimeter energy contained 
within a cone of radius R = ^(Ar/) 2 + (A</») 2 = 0.4 
around the electron direction excluding the energy asso- 
ciated with the electron. A more detailed description of 
the central electron selection can be found in hj. The 
forward electrons are selected by requiring an isolated 
energy cluster with E T > 20 GeV, the ratio of energy de- 
tected in the HAD and EM calorimeters to be less than 
0.05. The tracks are reconstructed using COT informa- 
tion in the region \rj\ < 1.6, while at higher \rj\ tracks are 
reconstructed using the SVX detectors alone. In order to 
reduce the charge misidentification and backgrounds, ad- 
ditional requirements for the forward tracks are imposed 
such as requiring the extrapolated charged-particle posi- 
tion to be consistent with the position measured in the 
calorimeter. Candidate W — > ev events are required to 



5 



have exactly one as well as i£ T > 25 GeV. The fi- 
nal W — •> ev data sample contains 537,858 events with 
a central electron and 176,941 forward electron events. 
To evaluate the detector acceptance and resolution for 
W — ► ev events we use the pythia [l8[ event generator 
followed by the CDF detector simulation. 

We determine the neutrino's longitudinal momentum, 
within a two-fold ambiguity, by constraining the ev mass 
to be that of the W boson. This ambiguity can be re- 
solved on a statistical basis from the known V — A decay 
distribution using the decay angle between the electron 
and the proton in the W rest frame, 8*, and from the 
W + and W~ production cross sections as a function of 
W rapidity {da ± /dyw)- To do this we assign a weight- 
ing factor to the two rapidity solutions, depending on the 
charge of the W boson, ii\ 2 '■ 

± J , ±(cosflf i2j y 1)2j pSS r )tr t (yi, 2 ) 
1,2 P± (cos 61 , yi , p*>± (j/i ) + P ± (cos 8* 2 ,y 2 , j#>± (y 2 ) 

(2) 



where 



P±{cos8*,y w ,p¥) = (lTco S 0*) 2 +Q{y w , p W)(l±co S 0*) 2 . 

(3) 

The ± signs indicate the W charge and the indices 1 
and 2 are for the two W rapidity solutions. The dif- 
ferential cross section as a function of yw is determined 
using a next-to-next-to-leading order (NNLO) QCD cal- 
culation using the MRST 2006 NNLO PDFs 
The ratio of the (1 + cos#*) 2 to (1 — cos6>*) 2 angular 
distributions, Q(yw ,Pr), m Eq. [3]is determined by the 
quark versus anti-quark composition of the proton using 
the event generator mc@nlo [2l[. This ratio is evalu- 
ated as a function of yw and the W transverse momen- 
tum . Although the weighting factor given by Eq. [5] 
depends primarily on the W + and W~ cross sections, it 
does have some weak dependence on the input W charge 
asymmetry. Therefore, this method requires us to iterate 
the procedure to eliminate this dependence. 

Correct charge identification is crucial for the measure- 
ment of the charge asymmetry measurement, because 
it directly affects the yield for a particular charge and 
yw and is corrected for on an event-by-event basis. The 
charge misidentification rate (Charge MisID) is measured 
as a function of rj using Z — > ee events where both elec- 
trons are identified as having the same charge sign. The 
misidentification rate ranges from (0.18 ± 0.05)% for 
\r)\ < 1.1 to (17.26 ± 2.02)% for > 2.04. 

The A{yw) values are corrected for the backgrounds to 
W — > ev candidates. We consider W — > tv, where the r 
decays leptonically to an electron plus neutrinos, as con- 
tributing to the signal and is included in the overall signal 
acceptance. The background fractions due to Z — > e + e~ 
events where one of the electrons is not reconstructed 
and mimics a neutrino are (0.59 ± 0.02)% for central 
electrons and (0.54 ± 0.03)% for forward electrons. The 



small contamination from the Z — > t + t~ process is found 
to be (0.10 ± 0.01)% for both central and forward elec- 
trons. The background from misidentified jets (QCD) 
is estimated by fitting the isolation distribution of elec- 
trons. Electrons in the calorimeter are characterized by 
having most of their energy deposited within an isola- 
tion cone centered on the electron, while jets may have 
significant energy deposits outside this cone. The QCD 
background fraction for central and forward electrons are 
(1.21 ± 0.21)% and (0.67 ± 0.18)%, respectively. 

The scale and resolution of the electromagnetic 
calorimeter energy and the missing transverse energy can 
affect the measured W rapidity and thus the asymme- 
try measurement. The EM calorimeter energy scale and 
resolution are tuned in the simulation to reproduce the 
Z — > e + e~ mass peak in data. The uncertainties on 
the energy scale and resolution for central electrons are 
measured to be ± 0.05% and ± 0.07%; for forward elec- 
trons they are ± 0.3% and ± 0.8%, respectively. The 
hadronic showering, the boson recoil-energy, and the un- 
derlying event energy of the hadronic calorimeter energy 
measurement play important roles in determining the $ T . 
The simulation for the calorimeter deposition in W — > ev 
events is tuned to provide the best possible match with 
data, including its dependence on r\. The uncertainty on 
the transverse recoil energy scale is ± 0.3% and ± 1.4% 
for central and forward electrons, respectively. 

We also investigate potential sources of a charge bias 
and r\ dependence in the kinematic and geometrical ac- 
ceptance of the event (estimated with simulated data) 
and efficiencies of the trigger and the electron identifica- 
tion (measured with data). The trigger efficiencies for the 
central and forward electrons are measured using data 
from independent triggers. We find the trigger efficien- 
cies do not depend on charge, but do depend on the rj 
and E T of the electron. The average trigger efficiencies 
for the central and forward electrons are (96.1 ± 1.0)% 
and (92.5 ± 0.3)%, respectively. Electron identification 
and track matching efficiencies (ID) are measured in data 
and simulation using the Z — > e + e~ channel. 

The choice of PDF sets has an effect on the shape of the 
da^/dyw distribution, as well as on the ratio of quarks 
and anti-quarks in the angular decay distribution. We 
use the 40 CTEQ6.1 error PDF sets [22j and re-determine 
the da /dyw production cross section and the angular 
distribution of cos 6* for each error PDF set. 

As expected, the data are found to be invariant un- 
der CP transformations A(yw) = — A{— yw), the two 
sets of points are in statistical agreement, so we com- 
bine the A(yw) bins with the complementary — A(— yw) 
bins in order to improve the precision of this measure- 
ment. We quote the statistical combination of A(yw) 
with — A{— yw), usin g th e Best Linear Unbiased Esti- 
mate (BLUE) method [23j , accounting for all correlations 
for both positive and negative yw bins. The statistical 
correlation coefficient between bins is found to be < 0.05. 
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TABLE I: Statistical and systematic uncertainties for the W production charge asymmetry. All values are (xlO 2 ) and show 
the correlated uncertainties for both positive and negative rapidities. 



\yw\ 


Charge 


Back- 


Energy Scale 


Recoil 


Electron 


Electron 


PDFs 


Stat. 




MisID 


grounds 


& Resolution 


Model 


Trigger 


ID 






0.0 - 0.2 


0.02 


0.04 


0.01 


0.11 


0.03 


0.02 


0.03 


0.31 


0.2 - 0.4 


0.01 


0.09 


0.04 


0.22 


0.08 


0.07 


0.08 


0.32 


0.4 - 0.6 


0.02 


0.11 


0.06 


0.22 


0.13 


0.17 


0.15 


0.33 


0.6 - 0.8 


0.03 


0.15 


0.07 


0.34 


0.14 


0.30 


0.22 


0.32 


0.8 - 1.0 


0.03 


0.20 


0.07 


0.42 


0.11 


0.47 


0.24 


0.34 


1.0 - 1.2 


0.04 


0.18 


0.08 


0.33 


0.09 


0.69 


0.27 


0.38 


1.2 - 1.4 


0.05 


0.18 


0.15 


0.67 


0.06 


0.78 


0.28 


0.43 


1.4 - 1.6 


0.04 


0.14 


0.14 


1.10 


0.04 


0.85 


0.28 


0.50 


1.6 - 1.8 


0.08 


0.12 


0.26 


0.92 


0.03 


0.89 


0.29 


0.55 


1.8 - 2.05 


0.22 


0.13 


0.31 


0.82 


0.06 


0.80 


0.34 


0.62 


2.05 - 2.3 


0.44 


0.21 


0.53 


0.59 


0.17 


0.85 


0.42 


0.83 


2.3 - 2.6 


0.45 


0.19 


0.62 


0.40 


0.27 


0.86 


0.50 


1.10 


2.6 - 3.0 


0.14 


0.10 


0.60 


0.43 


0.28 


0.65 


0.53 


2.30 



Table U summarizes the statistical and systematic uncer- 
tainties on 

The measured asymmetry which combines 

the positive and negative yw bins, is shown in Fig. [1] 
Also shown are the predictions of a NNLO QCD cal- 
culation using the MRST 2006 NNLO PDF sets and a 
NLO QCD calculation using the CTEQ6.1 NLO PDF 
sets, which arc in agreement with the measured asym- 
metry. Values of A(yw) and the total uncertainty for 
each \yw\ bin are listed in Table HT1 Since this measure- 
ment depends on the width of the W, in particular for 
the highest yw bin, the bin centers account for the W 
rapidity and W mass range accepted in each bin. We 
correct the bin centers to the value of < \yw\ > (average 
of W + and W~ rapidities) for which the asymmetry is 
equal to the one for a fixed W mass of 80.403 GeV/c 2 . 

In conclusion, using a new analysis technique we re- 
port the first direct measurement of the W boson charge 
asymmetry from Run II of the Tevatron, using data from 
1 fb -1 of integrated luminosity taken with the CDF II de- 
tector. Since the total uncertainties are smaller than the 
uncertainties coming from PDFs, as is also shown in ||, 
this direct measurement of the asymmetry is more sen- 
sitive to the ratio of d/u momentum distributions in the 
proton at high x than previous lcpton charge asymme- 
try measurements. This result is therefore expected to 
improve the precision of the global PDFs fits. 
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FIG. 1: The measured W production charge asymmetry 
and predictions from (a) NLO CTEQ6.1 and (b) NNLO 
MRTST2006, with their associated PDF uncertainties. 
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TABLE II: The W production charge asymmetry with total 
systematic and statistical uncertainties. 



\yw\ 


< \yw\ > 


A{yw) 


&sys 


O sys-\-stat 


0.0 - 


0.2 


0.10 


0.020 


±0.001 


±0.003 


0.2 - 


0.4 


0.30 


0.057 


±0.003 


±0.004 


0.4 - 


0.6 


0.50 


0.081 


±0.004 


±0.005 


0.6 - 


0.8 


0.70 


0.117 


±0.006 


±0.006 


o o 
0.8 - 


l.U 


0.89 


0.146 


±0.007 


±0.008 


1.0 - 


1.2 


1.09 


0.204 


±0.008 


±0.010 


1.2 - 


1.4 


1.29 


0.235 


±0.011 


±0.012 


1.4 - 


1.6 


1.49 


0.261 


±0.014 


±0.015 


1.6 - 


1.8 


1.69 


0.303 


±0.014 


±0.014 


1.8 - 


2.05 


1.91 


0.355 


±0.013 


±0.014 


2.05 


- 2.3 


2.15 


0.436 


±0.013 


±0.016 


2.3 - 


2.6 


2.40 


0.537 


±0.014 


±0.018 


2.6 - 


3.0 


2.63 


0.642 


±0.012 


±0.026 
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